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ABSTRACT 

VLA, MERLIN and Hubble Space Telescope imaging observations of the extended 
regions of the symbiotic system CH Cygni are analysed. These extensions are evidence 
of a strong coUimation mechanism, probably an accretion disk surrounding the hot 
component of the system. Over 16 years (between 1985 and 2001) the general trend is 
that these jets are seen to precess. Fitting a simple ballistic model of matter ejection 
to the geometry of the extended regions suggests a period of 6520 ± 150 days, with 
a precession cone opening angle of 35 ± 1 degrees. This period is of the same order 
as that proposed for the orbital period of the outer giant in the system, suggesting 
a possible link between the two. Anomalous knots in the emission, not explained by 
the simple model, are believed to be the result of older, slower moving ejecta, or 
possibly jet material that has become disrupted through sideways interaction with 
the surrounding medium. 
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1 INTRODUCTION 



system are likely to be coplanar and eclipsing, with eclipses 



Symbiotic stars occupy an extreme and relatively poorly 



understood region of the binary c lassification scheme. The 
name -vtas comed by LVierril (1941) to describe stars wiiicii 
appeared to have a combination spectrum: that of high exci- 
tation lines usually associated with a hot, ionised nebula su- 
p erimp hsed on a cool continuum witii prominent absorption 
features consistent with a late-type star. At present they are 
understood to be interacting binaries (with orbital periods of 
a few to tens of years) consisting of a cool giant losing mate- 
rial mostly via the stellar wind and a hot, luminous compact 
object which ionises a portion of the cool component wind 
(seaquist et al. 1984). Such a state of affairs represents the 
so-called quiescent phase, which can be interrupted by pe- 
riods of activity. The active phases start with an eruption 



of the liut stai, an event indicated pliuLuiiieLiically by an 
increas e of the atar'3 brightnG33 by 2 - 6 mag, and/or apcctro — 
scopically by high velocity broad emission features from the 
central star. Both radio and Hubble Space Telescope (HST) 
imagin g can directly resolve t he remnants of such dramatic 
events (|Hack fc Parcscc 199^ ^ E^res_et^aL^^995 ; Kenny et al 
199e| ; [Richards et al. 1999| ; [Eyres et al. 2001bp . 



separated by around 5700 days. Later studies (Mikolajewski 
et al. 1990b| and references therein) confirmed this period 



Further spectral studies of the system (Hinkle et al 



199J ) suggested that CH Cyg is probably a triple-star system 
consisting of an inner 756-day period binary which is orbited 



by an unseen G-K dwarf on a 5300 day orbit. Skopal et al 
(1996| ) uncovered 756 day interval eclipses which show that 
all three stars in the system are likely to be in coplanar 
orbits. 

Each of the three outbursts seen since 1978 was accom- 
panied by high velocity broad emission features consistent 
with mass outflows. D uring 1984-85 the material wa s ejected 
at -600-2500 km s"^ ( [Mikolajewski fc Tomov IQSej ), whilst 

Doradic and 



the 1992-95 active phase was characterized by s 
in part bipolar outfiow at ~1000-1600kms~^ (Leedjarv & 



Mikolajewski 199(; 3kopal et al. 1996 [ijima 1996) and fi- 



nally, during the recent, 1998-2000, out burst mass outfiows 
at about 1000 km s~^ were observed ( ^yres et al. 2001a , 
hereafter Paper II). The outfiows may be correlated with a 
significant increase of the radio emission and the radio light 
curves during thes e periods fit well with the optical ones 
( Kenny et al. 199(: ) . The 1984 mass ejec tion has been linked 
to the emergence of bipolar emission (Taylor et al. 198f) 
which has been attributed to high velocity (~1000kms~ ) 
ski et al. (1987 ) showed that the orbits of the stars within the jets. Non-thermal emission features have been discovered in 



The symbiotic star CH Cygni displays partic ularly com- 
plex behaviour . Optical spectroscopic studies by Mikolajew- 
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these jets from radio observations (Crocker et al. 2001, here- 
after Paper I), explained by shocked regions that arise when 
the high-speed ejecta interact with existing wind material. 

This paper aims to extend our understanding of the 
jets in this interesting system by developing theories and 
models to explain the change in their morphology in the 
years following their emergence. 

In this paper we present analysis of radio data from the 
VLA and MERLIN, along with HST observations. 



Table 1. Restoring beams used at different radio observing fre- 
quencies. 



Instrument 


Frequency (GHz) 


Beamsize (arcseconds) 


VLA 


1.5 


1.5 




4.8 


0.5 




8.4 


0.25 




15.0 


0.15 


MERLIN 


1.5 


0.15 




4.9 


0.05 



2 OBSERVATIONS 

CH Cygni has been observed on several occasions since 1985, 
in several different radio wavelengths and optical domains, 
further details of which are given in Paper I. All observations 
have been presented previously with the exception of the 
2000 November 10? VLA data which is shown here for the 
first time. 

Radio data re duction was c arried out using the AIPS 
software package (Greisen 1999). The flux density of the 
phase reference source was determined, and the data cali- 
brated to correct phase and amplitude instabilities. These 
corrections were applied to the observations of CH Cygni 
and a map produced by deconvolving the instrumental re- 
sponse (the so-caUed "dirty beam" ) from the map using the 
well-established CLEAN algorithm ( Hogbom 1974 ). 

Observations made on 1999 August 12 with the Hub- 
ble Space Telescope's Wide Field and Planetary Camera 2 
through the F502N filt er, corresponding t o the [O HI] line 
at AA4956 and 5007A) (iBiretta et al. 200l|), are used here as 



Table 2. Position angle of the central region of the nebula of 
CH Cygni 



these, with a diffraction limted resolution of around 0.05 arc- 
seconds, offer direct comparison with the spatial resolution 
of the radio data (see Table 0. Each observation was made 
up of two sub-exposures, allowing cosmic ray subtraction to 
be carried out (see Paper II and references therein). 



3 INITIAL DATA ANALYSIS 

A preliminary examination of radio maps of CH Cygni 
(taken with the VLA and MERLIN), as well as the HST 
optical images (using the F502N [OHI] filter) show that the 
position angle of the extended regions changes drastically 
over the 15 years between the earliest (1985) and most re- 
cent (2000) observations. 

A more quantitative analysis requires a measurement of 
the position angle of the central regions of the radio nebula. 
A problem inherent in this method is that the geometry 
of the restoring beam (in the case of the radio observations) 
may dominate over any natural geometry of the region. This 
arises as a result of the {u, v) coverage of the instrument, 
which depends upon the altitude of the source above the 
horizon. 

To avoid these p roblems, the ra dio maps were recreated 
using the methods of Hogbom (1974) but a circular restoring 
beam was imposed upon them, as laid out in Table ^. 

Two dimensional Gaussians were fitted to the bright 
peak in the centre of each map, and the position angle (an- 
ticlock wise through east from north) of the semi-major axis 



I ns t r urn G n t 


Julis-ii D&tG 


Year 


1/ (Ghz) 


Position Syngls { } 


VLA 


2446087 


1985 


14 


137.7 ±0.3 


VLA 


2446509 


1986 


1.6 


135.4 ± 1.2 


VLA 


2446509 




4.8 


133.9 ± 0.3 


VLA 


2446509 




14 


123.2 ± 1.6 


VLA 


2446509 




22 


131.1 ± 10.7 


VLA 


2446576 




1.6 


136.7 ± 1.0 


VLA 


2446576 




4.8 


132.2 ± 0.4 


VLA 


2446576 




14 


117.2 ±0.7 


VLA 


2446576 




22 


113.3 ±4.4 


VLA 


2447003 


1987 


1.5 


134.9 ± 2.0 


VLA 


2447003 




4.8 


110.5 ± 0.5 


VLA 


2447003 




15 


91.8 ±1.2 


VLA 


2447450 


1988 


1.6 


138.9 ± 5.6 


VLA 


2447450 




4.8 


109.8 ±4.5 


VLA 


2447505 




14 


111.7 ±20.6 


VLA 


2448488 


1991 


4.8 


112.0 ±9.7 


MERLIN 


2448827 


1992 


4.8 


121.0 ±9.9 


MERLIN 


2449911 


1995 


4.8 


125.5 ± 5.0 


VLA 


2449916 




1.6 


138.9 ±4.8 


VLA 


2449916 




8.4 


132.0 ±2.0 


VLA 


2449916 




14 


138.0 ± 16.8 


VLA 


2449949 




1.6 


139.5 ± 3.1 


VLA 


2449949 




4.8 


134.9 ± 1.4 


VLA 


2449949 




8.4 


126.7 ± 1.1 


MERLIN 


2451235 


1999 


4.8 


174.3 ±4.7 


MERLIN 


2451307 




1.7 


166.5 ± 7.2 


HST 


2451402 




[OHI] 


162.2 ± 0.1 


VLA 


2451447 




4.8 


164.1 ±0.7 


VLA 


2451447 




8.4 


169.9 ±0.8 


MERLIN 


2451481 




1.6 


166.5 ± 7.2 


MERLIN 


2451607 


2000 


4.8 


200.5 ± 5.2 


MERLIN 


2451696 




4.8 


185.4 ±8.0 


VLA 


2451858 




1.6 


177.2 ± 5.2 


VLA 


2451858 




4.8 


171.5 ± 1.4 


VLA 


2451858 




8.4 


192.2 ± 1.0 


VLA 


2451858 




14 


197.6 ± 10.3 


MERLIN 


2451961 




4.8 


198.6 ± 18.7 



the AIPS task 'IMFIT', and are derived from the ratio of 
the axes of the fitted ellipse. These measurements are shown 
in Tabic H 

Since the details of the alignment of the CH Cygni sys- 
tem are not known, it was assumed that the orbital plane 
of the system is bei ng viewed edge-on, as ec li pses are seen 



of the r 3sulting ellipse was taken as the position angle of the 
central material. The uncertainties are those produced by 



within the system (Mikolajewski et al. 1987; Skopal et al 



1996). The extension might be expected to precess about an 
axis perpendicular to the orbital plane if this is associated 
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Figure 1. Variation of the position angle of the central region of CH Cygni between 1985 and 2000, along with the fitted curve described 
in Section pi The fitted curve should be treated with caution as the data set does not appear to cover a full period of precession. 



with the coUimating mechanism. Assuming that the rate of 
precession is constant, the position angle, 9pA should follow 
a simple sinusoid so that at time t the position angle is given 
by 



9pA = 6'mid + Socos 2tt X 



(1) 



Such a curve was fitted, giving a precession of period P — 
8373 ± 800 days about a position angle on the sky of ^mid ~ 
146.7±3.4 degrees. The opening angle of the precession cone 
was 6*0 = 37.1 ±4.7°. 

This precession period is longer than the 5294± 117 day 
period quoted for the length of th e orbit of the outer giant of 
the system by Hinkle et al. (1993). The precession period we 
see could be a resonance of this, the 756 day orbital period, 
or both. However, the period we find should be treated with 
some caution, as our data set does not appear to cover a full 
period of precession. 



4 EVIDENCE OF PRECESSION WITHIN THE 
RADIO JETS 

There is a large deviation of some of the points from the 
fitted curve. A flaw inherent in the method of fitting gaus- 
sians is that it will naturally pick up the position angle of 
material at a certain distance from the core - that distance 
being determined by the resolution of the map. Hence what 
is being picked up is the position angle of material that was 
ejected some time before the observation was made. 



If the material in the extension is being ejected ballisti- 
cally from a precessing source then it would be expected to 
trace out a conical helix of a form se en in many precessing 
objects (Zaninetti & van Horn 1988). Expansion has been 
detected in the jet since the first radio jets were seen by 
Taylor et al. (198e| ). 

Precessing jets are frequently seen in astronomical ob- 
jects, one of the best studied being SS433. Its extended neb- 
ula exhibits a characteristic "S" profile that has been at- 



tributed to a precessing jet (Hjellming & Johnston 1981) 
with great success. Code derived from the work of R. E. 
Spencer (priv. comm.) was adapted for CH Cygni by re- 
moving the relativistic behaviour of the jets. This code gives 
the on-the-sky positions of discrete knots in the jet at a 
given time, as determined by a set of parameters. These 
patterns were compared to the distribution of the 40 bright- 
est CLEAN components in each map at a given frequency 
and the deviatio n minimised through use of the Downhill 
Simplex Method (|Nelder fc Mead 1965| ). 

The process was carried out using the set of 5 GHz maps 
taken by the VLA in the highest resolution "A" configura- 
tion, as these provided the most extensive coverage of the 
system over time. 



4.1 Simple ballistic model 

The simplest model of radio jet precession is one in which 
material is ejected continuously, at a constant velocity, from 
the poles of a precessing object. The locus described by the 



4 M. M. Crocker et al. 



Table 3. Fitted parameters for the simple ballistic model of jet 
ejection and precession. 



Parameter 


Value 


P 


6519 ± 153 days 


To 


JD2444038 ± 128 


^mid 


140° ± 1° 


00 


35° ± 1° 


VQ 


1263 ± 18 km s"! 


i 


88° ± 1° 



points along the resulting curved path at a time t is de- 
scribed (see Figure ^ in parametric form by 

(2) 
(3) 
(4) 



X — Vo{t — t) COS^o 

y = vo{t " t) sin6o sin((/!>) 
z = vo{t — t) sin^o cos{(j}) 



where x, y and z are the coordinates of a given knot of 
material, ejected from the central source (at the origin) at a 
time r with velocity vq. The opening angle of the precession 
cone is 6o, measured from the x axis, and the azimuthal 
angle of the jet (relative to the y axis) is $ = Qt where f2 
is the precession rate. 

Such a model has been used to fit cur ves to the high - 
velocity jets of the microquasar SS433 by Spencer (1984). 
For this model, the parameters were: P, the precession pe- 
riod of the jet (= ^); tq, a fixed time corresponding to a 
point when the precession angle of the jets was along the 
line of sight (the y axis in Figure |^; dmid the position angle 
(on the sky) of the axis of precession; 9o, the opening angle 
of precession cone; wo, the velocity of the material ejected in 
the jet and i, the inclination angle of the precession axis to 
the line of sight. 

The size of the initial simplex was set up so that it 
had p aram eters centred aroun d P = 5294 days (iHinklel 
et al. l|99^), To =JD2446275 ( Klikolajewski et al 1990aD , 
Omid = 140°, do = 30° (Section |), V = 1230 kms"^ (Paper 
I) and i = 90°. The process was allowed to run until the 
standard errors of the evaluated fits became smaller than 
0.01. The "best-fit" parameters from this run were then fed 
in as the starting values for a second run. This process was 
repeated so that four runs were carried out. The final set 
of parameters is shown in Table ^ Results very similarb to 
this were found if the starting period was set to 8373 days, 
as found in Section ^. The resulting ephemerides, resulting 
from material being ejected in discrete bullets every 50 days 
and following the model described above, are shown in Fig- 
ures I and ^. 

The model fits reasonably well to the geometry seen 
in the maps, especially in those resulting from observations 
closely following the emergence of the radio jet in 1986. In 
these images (JD2446509 and JD2446576) the jets are two- 
sided and well extended, making fitting more robust. 

In later epochs, when the jet is more one-sided, the 
curves still fit the morphology seen in the radio maps well. 
The fact that such a fit is achieved so long after the initial 
outburst (and following at least two other outbursts) demon- 
strates that the behaviour of the bipolar jets is not far from 
a simple ballistic model of the kind applied successfully to 
SS433 and other sources. 




Figure 2. Coordinate system used to describe a simple precessing 
jet. The hot component is at the origin, and the jet is represented 
by discrete packets of material with coordinates (x, y, z). 



4.2 Disrupted jet model 

The greatest deviations of the simple model from the ob- 
servations are seen in the 1991 map (which heis a low peak 
brightness and high background noise level) and in the mo- 
tion of the NW component in the 1987 and 1988 maps. If 
the component is the same object in all of the maps, then it 
is moving in a fashion entirely different to a simple, ballistic, 
bullet of material. By inspection, it appears to be deceler- 
ating, perhaps due to interaction with dense material. 

The synchrotron theory put forward in Paper I required 
ejecta to be impacting upon previously ejected stellar wind 
material, in order to cause shocks and local condensations 
in the magnetic field. This might be expected to cause some 
degree of deceleration and disruption. 

Direct ev idence of decelera tion is provided by optical 
spectroscopy (Skopal et al. 2001), which detects the velocity 
of material within a few hundred R© of the central stars. 
Terminal velocities of up to 2500 km s~^ are seen. Radio 
maps. 



are seen. 

which observe the motion of material much further 
out (from a few AU up to 100s of AU), never detect on- 
the-sky velocities of more than 1500 km s~^, assuming a 



distance to CH Cygni of 268 pc (Viotti et al. 1998). This is 
evidence of a large degree of deceleration between these two 
scales. 

The actual offsets of the apparently-slowed component 
from the central radio peak at the different observing epochs 
are shown in Table ^. 

A c omponent moving ballistically, as described in Sec- 
tion 4.1, would exhibit a constant position angle (Spa) and 
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Figure 3. Fitted ballistic material ejection model overlaid upon VLA 5GHZ radio maps (1986-1988). 
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Figure 4. Fitted ballistic material ejection model overlaid upon VLA 5GHZ radio maps (1991-2000). 
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Table 4. Offsets of the anomalous northwestern component of 
the radio jet of CH Cygni. 



Date 



Band 



Offset (mas) 
RA Dec 



0PA (°) 



R (mas) 



2446509 U 

2446576 U 

2447003 C 

2447450 C 



671 ± 15 904 ± 13 

546 ± 47 902 ± 20 

773 ± 75 1195 ± 29 

532 ± 75 1430 ± 49 



323.4 ±1.0 1130 ±10 

328.8 ±2.8 1050 ±30 

327.1 ± 1.8 1420 ±30 

339.6 ±3.2 1530 ±50 



a radial distance from the core (R) that increased linearly 
with time. If there is simple deceleration, 9pA would remain 
constant while R changed according to an environment- 
dependent model. 

Cornpared with the ^ 1 ^.HO km .s'^ velocity .seen in thp 



main p ortion of the lets, the north western emission is al- 
most stationary on the sky. It is possible that the material 
has indeed almost stopped moving, perhaps due to interac- 
tions with the surrounding environment. The jet may have 
become entirely disrupted through shock-interactions and 
so would no longer obey the simple ballistic model. 

For a jet that has a time-variable ejection direction, 
such as a precessing bipolar jet, shocks will form between 
the jet material and the surrounding medium not only in 
the main direction of jet propagation but also as a re sult of 
the tran sverse interactions between the two regions ( Dyson 
19870 

The transverse velocity of the jet becomes an important 
factor at large distances from the point of ejection. 

A colimated outflow will only survive until the trans- 
verse shock has crossed the di ameter of the jet. To simplify 
things, the approximation of Raga et al. (1993) is used so 
that 



/Swofisin^o 



(5) 



Most of these values can be estimated from work already 
carried out. The jet radius, rj, can be estimated by measur- 
ing the size of the anomalous region in the 1986 VLA U 
band map. The region was found to be approximately circu- 
lar with a mean diameter of 0.18 arcsec. Assuming a distance 
of 268 pc this corresponds to a linear size of 6.9 x 10^^ m. 
The initial jet ejection velocity was found in Paper I to be 
1.23 X 10^ m s"^ which is in agreement with that found by 
fitting to the the purely ballistic ejection model. The rate of 
precession can be found from the fitted period of 6519 days 
to be n = 1.1 X 10~*rad s~^. The opening angle, do was 
estimated previously to be 35° ± 1°. 

The main unknown is th e relative density parameter, f3. 
Following R,aga et al. (1993), /3 is assumed to be constant, 
so that the ambient medium is entirely homogeneous. This 
may not be the case, as a constant spherical outfiow from a 
star is expected to produce an ambient density that varies 
as pe oc r~^, and the presence of a jet would be expected 
to cause disruptions in the local region. However, without 
more knowledge of t his region, a homogeneous medi um is as- 
sumed for simplicity. Raga & Noriega-Crespo (1993) derived 
a ratio of Pe/pj ~ 0.05 from line ratios in the Herbig-Haro 
object HH 111 V. If this applies to CH Cygni, then /3 is 
approximately equal to 0.2. 

Applying these parameters to Equation |^ gives a value 



for the time taken for the disruptive shock to cross the jet as 
~ 7 X 10^ s. The distance the material has covered in this 
time (along the x axis in Figure |^ will be, for the parameters 
discussed above, Xd = 311 AU. The actual distance of the 
observed anomalies in the radio jet in 1986-1988 is indeed 
around 300 AU. 

To determine whether this disruption is occurring, di- 
rect imaging of the shocked region is necessary. This has 
been done in the case of the Herbig-Har o objects HH46/47 
/[Rcipurth fc Heathcote 199l| ) and HH34 ( [Reipurth fc Heath- ' 



cote 1992). Ha maps were subtracted from [S II] maps (at 



wavelengths 6717A and 673lA) to produce images in which 
the leading bow shock and trailing jet shock could be seen 
separately. The ratio of the intensity of emission in the Ha 
and [S II] was then used to derive the rel ative densities in 



the jet mater ial and surrounding medium ( Raga & Noriega- 



Crespo 199J). Unfortunately, whilst Ha emission maps of 



CH Cygni exist (Paper II), no maps in [S II] are available. 



4.3 Models with reflective symmetry 

If the geometry of the jets is a result of them being gravita- 
tionally attracted to the outer giant, they would be expected 
to exhibit a morphology in which the position angle of the 
material is a function of distance from the ejection source, 
and i n which there is reflective sym metry about the central 
core (Zaninetti fc van Horn 1988). Such a model was ap- 



plied successfully to 3C449 by |Lupton fc Gott (1982| ). Such 
a pattern does not fit the observed extensions in the CH 
Cygni system, implying that the precession is due to effects 
within the coUimating mechanism, rather than post-ejection 
gravitational infiuences. 



4.4 Ballistic model with time—dependent ejection 
velocity 

CH Cygni is observed to undergo irregular outbursts with 
associated variable mass outfiow rates. Such a state of affairs 
renders invalid the assumption that what is being observed 
is a simple, continuous jet with constant ejection velocity. 
A jet that results from truly irregular outflows is difficult to 
model analytically using the methods described within this 
work. 

Recent work, as described in Paper II, has discovered 
what may be two types of mass outflow from the system. 
One is optically thin, at a velocity measured spectroscopi- 
cally to be 1200 km s~^ and the other is optically thick and 
highly irregular. There is evidence that the increase in mass- 
transfer from the cool component to the compact component 
that occurs near to times of periastron leads to higher lu- 
minosity of the compact component and, as a consequence, 

greater radiatively driven mass outflow. 

Assuming a distance of 268 pc (Viotti et al. 1998), the 



movement between 1986 and 1988 shown in Table ^ repre- 
sents an on-the-sky velocity of only ~ 200 km s~^. Extrap- 
olating this motion backwards gives the time of ejection of 
this material, assuming ballistic motion, as around Decem- 
ber 1978. This time is near the beginning of the large out- 
burst that lasted between 1977 and 1984. Minor outflows 



have been observed spectrally d uring the outburst (Hack 
et al. 1988 ; Anderson et al. 1980 ) at velocities of a few hun- 
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dred km s^^. These may coincide with the ejection of this 
anomalous packet of material. 

Recent HST ob servations (from October 1999) made by 
Corradi et al. (2001) show the existence of loops of material 
in the [Oil] nebula at a distance of ~ 1.5 and ~ 5 arcsec. 
Echelle spectroscopy of these loops places their kinematical 
axis at around 100 km s~^. Allowing for projection effects 
(i.e. assuming they are tilted at 37° to the plane of the sky) 
their kinematical ages are approximately 15.2 years for the 
inner loop and 51 years for the outer one, giving the times 
of ejection of the loops as August 1984 and February 1949 
respectively. The 1984 event fits in well with the emer gence 
of the radio jets first observed by Taylor et al. (1986) and 
the model proposed in thi s section, but no records o f a 1949 
outburst could be found (Mikolajewski et al. 1990a). 



5 POSSIBLE CAUSES OF PRECESSION 

Precession of bipolar jets is a co mmon phenomenon in as- 
trophysics (Fukue & Yokoo 1986). The precession is usually 
attributed to motion of the colhmating mechanism of the 
jets. 

A simple colhmating structur e is a circu mstellar torus 



or disk, such as that postulated by Fukue (1983). Flickering 
in the optical bright ness of CH Cygn i has been observed 



on sev eral occasions (Wallerstein 1968; Mikolajewski et al 



1990b) and has been explained by the orbital motion of a 



bright point on a disk surrounding the compact component. 

There are several mechanisms by which such a disk 
could be induced to precess, and four of these scenarios are 
examined here. Only one of these is able to cause precession 
at the rate that is observed, given suitable estimates for the 
appropriate parameters. 



5.1 Precession of a tilted accretion disk 

The gravitational infiuence of the nearby cool component, 
or perhaps the possible third star in the system, could force 
precessi on in this disk, and hence the jets. Following the 
work of Merritt fc Petterson (1980 ), the angular speed of a 
torus or disk, with radius r, undergoing forced precession as 
a result of the infiuence of a companion is 



T 2 



mil Tih + rric 



where rrih is the mass of the hot component, rric is the mass 
of the cool component, a is their separation and a is the an- 
gle between the equatorial plane of the torus and the orbital 
plane o f the stars. E stimates for the parameters, from the 
work of ^kopal (1997| ), are mi, = O.2M0, nic = O.93M0 and 
a = 1.5 AU. This relation assumes that the disk is centred 
on the hot component, and the companion star causing the 
forced precession is the cool giant on the 756 day orbit. It 
would also be sensible to assume that a < 45° (Mikolajewski 



et al. 1990b) and < r < a. To obtain a precession period 



of 8373 days requires a disk radius of around 0.7 AU. This 
is a very large accretion disk, and it would be hard to imag- 
ine how such a disk could form purely via wind accretion 
mechanisms, making this scenario unlikely. 



If instead it is assumed that CH Cygni is a binary sys- 
tem on a 5294 day orbit, with a — 7.9AU the disk radius 
is found to be around 21 AU, around twice the size of the 
proposed semimajor axis of the outer giant's orbit. 



5.2 Slaved disk 

An attempt to explain the observed 164 da y period preces- 



sion of the b ipolar jets in SS433 was made by van den Heuvel 



et al. (198C ), by assuming that the compact component was 
surrounded by an accretion disk that precessed at a rate 
governed by the apsidal motion of the parent star. 

Using Equations 12, 13 and 16 of that paper, taking 
reasonable ste llar and orbital parameters of CH Cygni from 



Skopal (1997) a nd assuming that th e coefficient of apsidal 



motion is 0.005 (Schwarzschild 1958) suggests that the pre- 
cession period of a slaved disk would be around 12000 years. 
This result is a strong function of the radius of the cool com- 
ponent, so an increase in this star's size by a factor of 4 would 
reduce this precession period to about 15 years. 

To evaluate the possibility of such a model, precise 
eclipse timings need to be made in the future. 



5.3 Irradiated disk 

An accretion disk that is irradiated by a central source can 
become warped as a result of a torque th at arises due to 
radia tion that is re-emitted from the disk (Livio & Pringle 
199^). 



Such a warped disk might be expected to precess. In 
contrast to the gravitationally induced precession mentioned 
above, this is "self precession", i n that no companion s tar is 
required. Following the work of Livio & Pringle (1996), the 



period of this precession is of the same order as the time it 
takes for the warping to develop in the disk. Using sensible 
parameters for the mass of a possible accretion disk (around 
10~^ M0), this timescale is of the order of ~2000 years, 
making this an unsuitable model for precession. 



5.4 Magnetic precession 

The inner region of an accreti on disk surrounding a magn etic 
star, such as a white dwarf (Mikolajewski et al. 1990b), is 
subject to a twisting torque and a warping torque. These 
arise as a result of the interaction between the surface cur- 
rent of the disk and the horizontal m agnetic fi eld that is a 
feature of the dipole field of the star (Lai 1999). 



The warping torque arises through the twisting of the 
field that threads the disk vertically. The precessing torque 
is a result of the azimuthal screening current that occurs if 
the disk is diamagnetic. In this case, the vertical field lines 
cannot penetrate the disk. An azimuthal current is induced 
in the disk. The radial component of the central star's mag- 
netic field interacts with this current to set up a vertical 
force acting on the disk. This force is uneven across the 
disk, so a net torque results. This torque acts to cause the 
disk (and hence the coUimation mechanism of the jets) to 
precess about the magnetic axis of the star. 

Lai (1999) gives the angular precession rate of such a 
disk in Equation 2.35 of that paper. Presently, most of the 
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required quantities are unknown, so only an order of mag- 
nitude estimate is possible. The magnetic moment of the 



white cwarf was estimated to be 2 x 10 T m using the 



This would indicate the presence of a minor ejectio n event 
prior to the main ejection of the jets first seen by Taylor 



method of postno^7l98^ ) ; tli^ surface dens ity of the 
disk was assumed to be of the order 10 kg (Lin et al 



et al. (1986). The variable ejection vel ocity model discussed 
here and the Echelle spectroscopy of Corradi et al. (2001) 



198^)~]rhe angle between the disk axis and stellar rotation 
axis was taken to be equal to the opening angle of the pre- 
cession cone, 35°, found previously. The magnetic axis of the 
system was assumed to lie along the disk axis. 

The magnetically induced precession period is then di- 
rectly related to the radial distance, r, of the warped torus 
from the white dwarf. Since this varies as r^, small changes 
in r greatly affect the period. A warped torus at around 2 Rq 
from the star will precess with a period of around 6000 days. 

A magnetic coUimation and ejection mechanism for the 
j ets gives a possible exp lanation for their one-sidedness 
( Chagelishvili et al. 19961). Jets powered by a central mag- 



netic rotator QMikolajewski et al. 1990h| ) can become inher- 
ently asymmetric, if the solution to the magnetohydrody- 
namic equations for each hemisphere of the star are asym- 
metric. Such conditions arise when there is no change in 
the relative sign of the poloidal and toroidal fields across 
the equator of the star. The different magnetic states above 
and below the equator can then easily lead to a condition 
in which only one supports the coU imation and ejection of 
material via a jet (Wang et al. 199S). 



6 CONCLUSION 

It is certain that the bipolar jets of CH Cygni exhibit pre- 
cession, with a period controlled by the motion of the mech- 
anism resposible for the coUimation of the outfiow. In ad- 
dition, the variability of the mass-loss rate from the cool 
component causes the material to be ejected at a rate that 
is highly variable (between ~ 500 and 2000km s~^), leading 
to complications in predicting the outfiow geometry. The ac- 
tivity of the bipolar ejection seems to be tied to the state 
of the hot-component, with fast jets seen at or soon after 
times of optical outburst. This activity then decreases dur- 
ing times of quiescence. 

Although the precession period found here is similar 
to the orbital period of the outer giant in the system, the 
current time resolution of the observations does not rule out 
shorter periods. A precession that resulted from the motion 
of the inner, symbiotic, pair and had a similar period to its 
756 day orbit would need observations to be taken at least 
every year for it to be detected. 

The simple ballistic model fits the geometry of the neb- 
ula extremely well. Knots seen in the larger-scale radio maps 
(the VLA in C, X and U band) that do not fit the simple 
ballistic model can be explained as regions of jet disruption 
caused by sideways motion of the ejected material brought 
about as a result of the precession. The predicted distance 
at which this disruption would occur is in close agreement 
with the observations, although it relies upon several as- 
sumptions about the nature of the ambient material. Con- 
firmation of this model would require simultaneous optical 
images in both Ha and [S II]. 

An alternative explanation is that the anomalous mate- 
rial was ejected at a much slower velocity and at an earlier 
time than the majority of the gas in the bipolar nebula. 



suggest the emission of material with a velocity of no more 
than a few 100s of km~^ just before the main ejection of the 
high-speed material. 

The cause of the precession is unknown but, given re- 
alistic estimates for the masses and separation of the stellar 
components, warping of the colUmating accretion disk by a 
magnetic white dwarf is the only mechanism that can give 
precession periods similar to those found by model fitting. 
The other possible causes give periods that are several or- 
ders of magnitude too large. The magnetic explanation also 
allows for the one-sided nature of the emission seen in all 
VLA radio maps following the initial outburst in 1984. 
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